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^ . We apply a rigorous statistical analysis to the Constrained MSSM to derive the most probable 

ranges of the diffuse gamma radiation flux from the direction of the Galactic center and of the 
positron flux from the Galactic halo due to neutralino dark matter annihilation, for several 
I different choices of the halo model and propagation model parameters. We find that, for a 

. specified halo profile, and assuming flat priors, the 68% probability range of the integrated 7- 

\ ray flux spans about one order of magnitude, while the 95% probability range can be much larger 

\^ . and extend over four orders of magnitude (even exceeding five for a tiny region at small neutralino 

\ mass). The detectability of the signal by GLAST depending primarily on the cuspiness of the 

. halo profile. The positron flux, on the other hand, appears to be too small to be detectable by 

' PAMELA, unless the boost factor is at least of order ten and/or the halo profile is extremely 

cuspy. We also briefly discuss the sensitivity of our results to the choice of priors. 
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1 Introduction 



There is currently much evidence for the existence of large amounts of dark matter (DM) in the 
Universe. While its nature remains unknown, DM is likely to be made up of an exotic species 
of weakly interacting massive particles (WIMPs). A particularly popular WIMP candidate is the 
lightest neutralino x of effective low-energy supersymmetry (SUSY), which is stable due to R- 
parity [HE]- In addition to collider searches for SUSY and direct detection (DD) searches for a 
cosmic WIMP, a promising strategy is that of indirect detection (ID), i.e., a search for traces of 
WIMP pair-annihilation in the Milky Way. Since the annihilation rate is proportional to the square 
of the WIMP number density, of particular interest are the Galactic center (GC) and nearby clumps 
in the halo where the density of DM is believed to be enhanced. The aim of this paper is to provide, 
for the first time, a statistical measure for the prediction of 7-ray and positron signatures in low- 
energy SUSY, thus allowing one to assess high-probability regions for DM-annihilation signatures 
that could be observed by the GLAST (in orbit since June 2008) and PAMELA (launched 2006) 
satellites. Existing data from EGRET suggest a spectrally distinct excess of 7-rays up to ~ 10 GeV 
and the HEAT data indicate a possible excess in positron flux between 5 to ~ 30 GeV. GLAST 
and PAMELA will provide an order of magnitude more sensitivity. 

In assessing detection prospects of WIMPs there are two main sources of uncertainties. One 
comes from the underlying particle physics model where WIMP mass and annihilation cross section 
can vary over a few orders of magnitude. The other is astrophysical in nature and stems from 
substantial uncertainties in the DM distribution, both locally (local DM density and the existence of 
clumps) and towards the GC. Since the general Minimal Supersymmetric Standard Model (MSSM) 
suffers from a lack of predictability due to a large number of free parameters, it is interesting and 
worthwhile to assess WIMP detection prospects in more constrained and more well-motivated low- 
energy SUSY models, among which particularly popular is the Constrained MSSM (CMSSM) [3], 
which includes the minimal supergravity model [3]. By applying a statistical approach, we derive in 
the CMSSM most probable ranges of fluxes, thus bringing under control all the uncertainties of the 
particle physics side of WIMP detection. This is a major improvement over existing methods which 
are usually limited to the consideration of a few representative choices of points or slices in the 
parameter space. Detection prospects then become a function of specific astrophysical uncertainties 
only. 

In this Letter we employ a Bayesian Markov Chain Monte Carlo (MCMC) technique to efficiently 
explore the multi-dimensional parameter space of the CMSSM, and to include all relevant sources of 
uncertainty on the particle physics side [5l |6] (for a similar study, see [7] ) . Our Bayesian approach 
allows us to produce probability maps for all relevant observable quantities, thus establishing a 
complete set of predictions of the CMSSM. 

2 Bayesian analysis of the CMSSM 

The CMSSM is described in terms of four free parameters: a ratio of Higgs vacuum expectation 
values tan/3, and common soft SUSY-breaking mass parameters of gauginos, mi/2) scalars, mo, 
and tri-linear couplings, Aq. The parameters "io and Aq are specified at the GUT scale, 

Mgut — 2x 10^^ GeV, which serves as a starting point for evolving the MSSM renormalization group 
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equations for couplings and masses down to a low energy scale MgusY = ^m^_^ m^^ (where mj^ 
denote the masses of the scalar partners of the top quark) , chosen so as to minimize higher order loop 
corrections. At MsusY the (1-loop corrected) conditions of electroweak symmetry breaking (EWSB) 
are imposed. The sign of the Higgs/higgsino mass parameter however, remains undetermined. 
Here we set /i > 0. 

In deriving predictions for the observable quantities, one also needs to take into account the 
uncertainty coming from our imperfect knowledge of the values of some relevant Standard Model 
(SM) parameters, namely the pole top quark mass, Mt, the bottom quark mass at rrib, mh{mi))^^ , 
and the electromagnetic and the strong coupling constants at the Z pole mass Mz, aem{Mz)^^^ 
and as{Mz)^^^ , respectively (the last three quantities are all computed in the MS scheme). These 
four "nuisance parameters" are the most relevant ones for accurately predicting the SUSY spectrum 
and its observable signature. In our analysis we thus consider an 8-dimensional parameter space 
spanned by the above four SM and the four CMSSM parameters. 

In general, the results of a Bayesian analysis are expressed in terms of a posterior probability 
distribution (or more briefly, "a posterior"). By virtue of Bayes' theorem, the posterior is the 
product of the prior and the likelihood. The prior expresses the state of knowledge about the 
parameters before seeing the data, while the likelihood encodes the information coming from the 
observations (for further details on the Bayesian framework, see e.g. [9]). If the constraining power 
of the data is strong enough, then the posterior is effectively dominated by the likelihood and the 
prior distribution becomes irrelevant. However, if the observations are not sufficiently constraining, 
the posterior will retain a prior dependence. Therefore it is important to check to which extent the 
results based on the posterior pdf show a prior dependency. There are reasons to believe that for 
the CMSSM present data are not sufficiently powerful to completely override the prior, see [TO] for 
a detailed study of this issue. 

In our analysis we assume flat priors on both SM and CMSSM parameters over wide ranges 
of their values, encompassing the focus point region [8]. However, below we will comment on how 
our result change when one employs a flat prior in logiQmi/2 and logj^g'^o instead (which we call 
in the following "the log prior" for brevity). The reason for this alternative choice of prior is that 
that they are distinctively different. In particular, the log prior gives equal a priori weights to all 
decades for the parameters. So the log prior expands the low-mass region and allows a much more 
refined scan in the parameter space region where finely tuned points can give a good fit to the data 
(see |10) for details). Other choices of priors are possible and indeed physically motivated, and will 
be considered in future work. A recent discussion of some alternative prior choices in the CMSSM 
(motivated by considerations of naturalness and fine tuning) can be found in Ref. [11]. 

At every point in parameter space, we compute a number of observable quantities, and compare 
their values with the observational data listed in Ref. [GjH where also a detailed description of our 
procedure is given. We include all relevant collider limits, including direct limits on Higgs and 
superpartner masses, rare processes BR[B ^s7), BR[Bs /i"*"//") and recently measured Bg 
mixing, AMb^ , electroweak precision data {mw and sin^ ^efr) and the relic abundance of the lightest 
neutralino ^i-^/i^ assumed to be the cold DM in the Universe. We then use our MCMC algorithm |12j 

*We employ the WMAP 3-year relic abundance value assuming that neutralinos are the only dark matter compo- 
nent. Using the WMAP 5-year value instead would not change visibly our results. 
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to produce, for a given model of DM distribution in the Galactic halo, probability distribution maps 
in parameter space and various observables, including ID ones which are computed with the help of 
DarkSusy [13]. As we have emphasized in Ref. [6], current constraints, especially from 6 — > 57 favor 
the focus point region of large rriQ > 1 TeV and not so large mi/2 ^ 1-5 TeV (with mi^2 ^ 2.5mo). 



3 Gamma-ray flux from the Galactic center 

The differential diffuse 7-ray flux arriving from a direction at an angle ip from the GC is given by [2] 



where aiV is a product of the WIMP pair-annihilation cross section into a final state i times the 
pair's relative velocity and dN^/dE^ is the differential 7~ray spectrum (including a branching ratio 
into photons) following from the state i. Here we consider contributions from the continuum (as 
opposed to photon lines coming from one loop direct neutralino annihilation into 77 and jZ), 
resulting from cascade decays of all kinematically allowed final state SM fermions and combinations 
of gauge and Higgs bosons. The integral is taken along the line of sight (l.o.s.) from the detector. It 
is convenient to separate factors depending on particle physics and on halo properties by introducing 
the dimensionless quantity J^tp) = (l/8.5kpc) (0.3 GeV/cm^) dl p^{r{l,^)) [H]. The flux is 

further averaged over the solid angle AO representing the acceptance angle of the detector, and one 
defines the quantity J(AO) = (1/AO) /^^ J{^)dVt. 

Clearly, one of the crucial ingredients is the radial dependence of the WIMP density Pxi^)- Some 
popular profiles can be parameterized by [2] 

Px(0 = Po ^"'"'^ V. [1 + (^o/a)"]^ , (2) 
[l + (r/a)"] — 

where the halo WIMP density has been normalized to its local value, assumed to be /?o = 0.3 GeV/cm^. 
Table [T] gives the values of the parameters: a, a, /?, 7 and tq for some common choices. Here we 
consider the line-of-sight (l.o.s.) integration factor J in the direction of the GC, i.e., for tp = {). 
In the case of the cuspy profiles, in order to avoid a divergent behavior, we set a cutoff radius of 
Tc = 10~^ kpc. The total 7-ray flux from the cone AO centered on ifj and integrated over photon 
energy from an energy threshold i?th) is then given by 



$^(Afi)=/ dE^d<S>^/dE^{E^,An). (3) 



-5, 



With the launch of GLAST, which has angular resolution AO ~ 10 sr and sensitivity to fluxes 
larger than about 2 x 10~ cm~ s~ for photon energies E^ > 10 GeV [2T]|£| it is timely to investigate 
the global predictions of the CMSSM for a range of halo models. 



^Resolution and sensitivity in the range 30 MeV < < 10 GeV are energy-dependent and would require a more 
careful analysis. 
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Halo model 


a ( kpc) 


a /3 7 


J(lO-^sr) 


J(lO-'^sr) 


isothermal cored 


3.5 


2 2 


30.35 


30.40 


NFW 


20.0 


1 3 1 


1.21 X 10^ 


1.26 X lO'* 


NFW+ac 


20.0 


0.8 2.7 1.45 


1.25 X 10^ 


1.02 X 10'^ 


Moore 


28.0 


1.5 3 1.5 


1.05 X 10^ 


9.68 X 10^ 


Moore+ac 


28.0 


0.8 2.7 1.65 


1.59 X 10*^ 


3.12 X W 



Table 1: Parameters for some popular halo profiles: a spherically symmetric modified isothermal 
model [15], the Navarro, Frenk and White (NFW) model and the Moore, et al., (Moore) 
model |l7j. Everywhere tq = 8.0 kpc except for the isothermal case, where ro = 8.5 kpc. In the 
NFW and Moore, et al., models the effect of adiabatic compression due to baryons (marked with 
an additional +ac), is included. See also Ref. ^8]. 

Fig.[T]shows the joint probability distribution for the total flux from the GC above a threshold 
energy of 10 GeV versus m^, integrated over AJ7 = 10~^sr. The spread of values reflects the 
marginalization over all the four CMSSM and four SM parameters, thus fully accounting for all 
substantial sources of uncertainty on the particle physics side. Firstly, as can be seen from the 
figure, the 2-dim joint 68% probability range of lies between about 80 GeV and about 600 GeV. 
Secondly, for a given halo profile, and assuming flat priors in CMSSM parameters, we find that the 
68% probability range of is confined to lie within about one order of magnitude. On the other 
hand, the spread of the 95% probability range is much larger and at lower can extend over four 
or even five orders of magnitude. 

In order to examine the low mass region in more detail, we have redone our analysis for the log 
prior choice introduced above. As mi 12 and rriQ are the primary CMSSM parameters determining 
mass spectra of the neutralino, the other superpartners and the Higgs bosons, the log prior allows 
one to examine the low mass region in more detail, in particular by "expanding" the volume of the 
region 100 GeV < mi/2)"io 1 TeV. As we discuss below, the flat prior appears to produce an 
optimistic scenario as far as indirect detection signatures are concerned, while the log prior can give 
lower values of the fluxes and hence it leads to more pessimistic prospects for indirect detection. 
Ways of mediating between the two scenarios and to assess their relative plausibility will be explored 
in future work. 

Since the log prior gives more "weight" to lower values of both a-nd mo, not surprisingly, 
we have found that it leads to a large widening of mostly the lower boundary of the 68% probability 
range at low m^, while not affecting the flux ranges at larger values of the neutralino mass. For 
example, the 68% probability range widens to nearly three decades and, in the case of the Moore 
profile with adiabatic compression, can be as low as 1.2 x 10^^'^cm~^s~^ at ~ 100 GeV, but 
then it quickly raises and for > 200 GeV is not very different from the case of the flat prior. 
A more detailed discussion of the implications for CMSSM parameters of employing a log prior is 
given in Ref. [lOj . 

For a given prior, choosing a different halo profile merely amounts to shifting the total flux 
by the ratio of the values of J given in Table [H As expected, more cuspy profiles lead to higher 
predicted fluxes. We find that, in the CMSSM in the case of the Moore profile (with and without 
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Figure 1: The joint probability distribution for the 7-ray flux <I>^ from the Galactic center vs the 
neutralino mass for some popular halo profile models, assuming flat priors, as explained in the 



text. The dark/red (light/yellow) region shows the predicted 68% (95%) probability ranges for the 
Moore profile with adiabatic compression (for other profiles we indicate only the limits of the 68% 
region). We also plot the expected 5-o" detection threshold (neglecting background) for energies 
above E^y = 10 GeV for GLAST after 1 year of operation [21]. Notice that for GLAST to be able 
to detect the annihilation flux over the background one might require much larger fluxes than the 
sensitivity level plotted here 



adiabatic compression) and the NFW proflle with adiabatic compression, the continuum flux signal 
will be within the reach of GLAST, while for proflles with J(10~^sr) < 10^ it will not be detectable 
by GLAST. (We have checked that the case of ;U < and flat priors gives qualitatively similar 
results.) 

The differential 7-ray flux from DM annihilations is expected to exhibit a sharp drop-off in 
the energy spectrum as Er^ approaches m^. In Fig. [2] we plot 68% and 95% probability regions 
for the 7-ray differential flux for the NFW proflle, averaged over a solid angle AJl = 10~^sr (to 
allow a comparison with EGRET data), for the flat prior choice. Clearly, the current uncertainty 
on CMSSM parameters and hence on m-^ introduces a considerable spread in the predicted spectral 
shape of the signal. Additional uncertainty comes from the dependence on the priors. For example, 
for the log prior given above the 68% probability range of the differential photon range extends 
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Figure 2: Predicted 7-ray differential energy flux averaged over a solid angle Ail = 10^'^sr and 
fully accounting for current uncertainty in the CMSSM parameters, assuming flat priors. The 68% 
and 95% regions are for the NFW profile, all other cases can be obtained by rescaling them by 
the factors J given in Table [H Predictions appropriate for GLAST resolution (AO = 10~^sr) are 
obtained by dividing by 1.21 x 10^ and multiplying by the desired value of J(lO^^sr). We plot 
the expected GLAST 80" detection threshold (horizontal black/dashed line }21j). The three blue, 
dashed curves show sample spectra (for the values of specified in the figure) from our statistical 
scan. For comparison we plot EGRET diffuse data towards the GC (green squares [19]), EGRET's 
point-source subtracted fiux (blue empty squares ^20j) and H.E.S.S. (2004) data [23j (red squares) 
with 2(7 error bars. 



between 2.4 x 10~^^ GeVcm~^s~^sr~^ and 6.7 x 10~^ GeVcm~-^s~^sr~^. Thus, even if a positive 
signal were detected by GLAST, it would be difficult to infer from it the mass of the WIMP, 
especially at its lower values below some 200 GeV, with any reasonable accuracy. 



4 Positron flux from the Galactic halo 

Positrons can be produced either in direct DM annihilation, or from decays and hadronization 
of other products (gauge and Higgs bosons, etc), with the continuum spectrum from the latter 
usually dominating. Once produced, they propagate through the Galactic medium and their spec- 
trum is distorted due to synchrotron radiation and inverse Compton scattering at large energies. 
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bremsstrahlung and ionization at lower energies. The effects of positron propagation are com- 
puted following a standard procedure described in PH [22] , by solving numerically the diffusion- loss 
equation for the number density of positrons per unit energy dn^+/d£. The diffusion coefficient is 
parameterized as K{e) = i^o(3" + with Kq = 5.8 x lO^^ cm\-^ , a = 0.6 and e = E^+/l GeV, 
mimicking re-acceleration effects. The energy loss rate is given by 6(e) = te£^, with te = 10~^^s~^, 
and we describe the diffusion zone (i.e., the Galaxy) as an infinite slab of height L = 4kpc, with 
free escape boundary conditions. Changes in the above positron propagation model, especially 
K{£) (see e.g. [26l ES]), can potentially lead to variations by a factor of 5 to 10 in the spectral 
shape at low positron energy, < 20 GeV [27j. In this energy region the flux dependence on 
the halo profile is also substantial and, for the models in Table [H the flux can change by up to a 
few orders of magnitude (compare blue/dashed lines in Fig. [3]), since positrons from further away 
loose energy due to propagation. Most high-energy positrons, on the other hand, originate from the 
local neighbourhood the size of a few kpcs [25^ I29j. and their flux is less dependent of the halo and 
propagation dynamics. The flux can, however, be considerably enhanced by the presence of local 
DM clumps that survive merging processes and tidal stripping [28] , an effect that is usually param- 
eterized by a boost factor (BF), which can be of order 10. Recent studies have begun investigating 
the dumpiness dependence of the spectrum in more detail [30[ I29j. Finally, in order to reduce the 
impact of solar winds and magnetosphere effects on the model's predictions, it is useful to consider 
the positron fraction, defined as $e+/(*^e+ + '^'e-); where <I>e+ is the positron differential flux from 
WIMP annihilation, while <!>£- is the background electron flux. For background e~ and fluxes 
we follow the parametrization adopted in Ref. [25] from Ref. [26j. 

In Fig. [3] we show the predicted positron flux fraction in the CMSSM (for flat priors) for the 
NFW profile and a boost factor BF=1, alongside a compilation of observations, most notably from 
HEAT. Again, the uncertainty in the spectral shape is one of the main results of our analysis, which 
accounts for the current uncertainties regarding the CMSSM parameters. The 95% probability 
region peaks in the range 1 GeV < Ef,+ < 10 GeV, roughly in the region of the apparent HEAT 
positron excess, but the strength of the signal is insufficient for it to be detectable by PAMELA 
in the absence of a large boost factor [3T]. On the other hand, in that energy range the signal 
would be enhanced by more than one (two) order(s) of magnitude for a more cuspy profile such as 
the NFW (Moore) profile with adiabatic compression, as indicated in Fig. [3] for a sample spectrum 
corresponding to = 229 GeV. This is because in the case of more cuspy profiles more high- 
energy electrons coming from the GC are scattered to lower energies. For the case of the log prior 
on mi/2 and rriQ we again find a significant decrease of the lower boundary of the 68% range. For 
example, for the NFW profile with BF=1 the ratio in Fig. [3] can be as low as 3.8 x 10~^. 

We conclude that, for not exceedingly cuspy halo models, PAMELA is unlikely to be sensitive 
to positron fluxes in the CMSSM, since for the NFW profile the signal is more than two order of 
magnitude smaller than the background. This result would qualitatively hold even when taking 
into account the considerable uncertainties coming from the boost factor due to local clumps and 
changes in the positrons propagation model, each of which can potentially change the spectrum by 
up to a factor of 10. 
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Figure 3: Predicted positron flux fraction in the CMSSM. The 68% (dark/red) and 95% 
(Hght/yeflow) regions are for an NFW profile with a boost factor BF=1 and a specific choice of 
propagation model. We also show for comparison some of the current data. To illustrate the depen- 
dency of the spectral shape at low energies on the halo model, we plot the spectrum for the same 
choice of CMSSM parameters (with = 229 GeV) for three different halo models as indicated. 
In absence of a large boost factor, the signal appears too small to be detected by PAMELA. 

5 Summary 

In the framework of the Constrained MSSM, we have performed a Bayesian analysis of prospects 
for indirect dark matter detection via a diffuse 7-ray signal or a positron flux from the Galactic 
center. This has allowed us to provide a statistically rigorous assessment of the uncertainty from 
the particle physics side of the problem. 

We found that the prospects for GLAST to detect a diffuse 7-ray signal from the Galactic center 
depend primarily on the cuspiness of the DM proflle at small radii. For the choice of flat priors 
in the CMSSM parameters, the NFW model appears to be a borderline case, while a more cuspy 
halo would guarantee a signal for a 68% range of the CMSSM parameter space, except near the 
bottom end of the neutralino mass around 100 GeV, below the 68% probability range of m^. In the 
low mass region the sensitivity to the choice of priors remains however substantial. Adopting a log 
prior on ^i^d ^0 leads to a significant decreasing of the lower boundary of the 68% probability 
range of the 7-ray flux towards lower values at low niy. < 200 GeV, but at larger gives similar 
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results as with flat priors. 

On the other hand, a positron flux is unlikely to be detectable by PAMELA for both choices of 
priors, unless it is strongly enhanced by a nearby clump with a boost factor of at least of order ten. 
The latter conclusion is valid for a specific (although well motivated) choice of propagation model 
parameters. Assumptions regarding propagation parameters could however be easily relaxed in our 
framework. It would be straightforward to extend our treatment to include propagation model 
parameters as nuisance parameters and marginalize over them, as well. It is expected that such a 
procedure would increase the present, very substantial uncertainty as to the spectral shape, which 
we have shown is a consequence of the current lack of knowledge as to the preferred regions of the 
CMSSM parameters. Finally, it would also be interesting to repeat this analysis in a more general 
phenomenological SUSY model than the Constrained MSSM. While a richer phenomenology might 
help in explaining future signals should they be detected, it is also clear that a larger number of 
free parameters on the particle physics side will add to the difficulty of reliably predicting the shape 
and strength of both the 7-ray and the positron spectra. 
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